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An exact incompressible quantum liquid is constructed at the filling factor 1/m 2 in square lattice. 
It supports deconfined fractionally charged excitation. At the filling factor 1/m 2 , the excitation has 
fractional charge ±e/m 2 , where e is the electric charge. This model can be easily generalized to the 
integer lattice in any D dimensions, where the charge of excitations becomes ±e/m D . 



Quantum number fractionalization is one of the most 
intriguing phenomena in physics. The most astonish- 
ing discovery of the fractional excitation is arguably the 
Laughlin quasiparticle in the fractional quantum Hall 
effect (FQHE). Laughlin quasiparticle has charge e/q, 
where e is the electric charge and q is taken to be any 
odd integer. Different from quarks, it is indisputably a 
fractional excitation because 1) it is in the deconfined 
phase and 2) it is a new degree of freedom in addition to 
the ground state which is an electronic quantum liquid. 

In one dimension, quantum number fractionalization is 
a robust effect. There are many models in one dimension 
exhibiting the fractionalization of charge or spin in the 
excitation spectrum. On the contrary, it is difficult to be 
found in theoretical models and in experiments in the di- 
mensions greater than one, although there is no theorem 
to restrict its existence in higher dimensions. Whether 
or not deconfined fractionalization can in principle occur 
in any dimensions is still an open question. The under- 
standing to this problem is still primitive, neither is there 
a general theoretical construction. In this Letter, we an- 
swer the question positively by providing, to our best 
knowledge, the first microscopic model to host its exis- 
tence in the integer lattice in any dimensions. Specifi- 
cally, the two-dimensional case in the square lattice will 
be worked out in detail. The higher dimensional general- 
ization can be systematically constructed. It is the first 
theoretical model to host the exact ground state of in- 
compressible quantum liquid in any dimensional integer 
lattice. 

The model we consider here was originally proposed 
by one of us (Chern) and Lee pQ, where a lattice pair- 
hopping model of incompressible quantum liquid is con- 
structed in the (n — l)-dimensional triangular lattice. It 
is basically the generalization of the fractional quantum 
Hall effect in CP™ -1 , which is the 2(n — l)-dimensional 
complex projective space. In ref.p], we constructed the 
generalized Laughlin wavefunction in CP™ -1 and found 
the Hamiltonian to host it as the unique ground state 
which is well separated from the excitations by a finite 
energy gap. The "lattice" is the weight space of the co- 
herent states of SU(n), the isometry group of CP" -1 , 
in the (jp, 0, ...,0) representation, where p is an integer. 
Without losing generality, in the following we start with 



the simplest generalized case, namely CP 2 . After briefly 
reviewing the FQHE in CP 2 and taking the flat space 
limit from CP 2 to R 4 , we shall illustrate that the in- 
compressible quantum liquid can be constructed in the 
square lattice with periodic boundary. Furthermore, we 
shall show that the fractionally charged excitation, ap- 
pearing as a dislocation, is in the deconfined phase. Most 
importantly, our theory can be generalized to higher di- 
mensions easily. 

The Laughlin wavefunction in CP 2 can be written by 
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where u, v, w form the SU(3) fundamental triplet, N is 
the number of particles, m is any odd integer, and to- 
gether with p they label the SU(3) (mp, 0) representa- 
tion of which the weight space spans a two-dimensional 
triangular lattice with open boundary. The relation 
between N and p is given by N = d(jp, 0) , where 
d(p,0) = (p + l)(p + 2)/2 the dimension of the SU(3) 
(p, 0) representation. For general m, the filling factor 
v = d(p,0)/d(mp,0), which equals to 1/m 2 in the ther- 
modynamic limit. 

For convenience, we consider m — 3 in the follow- 
ing. It will be soon clear that it is not hard to gener- 
alize to any odd m. In the flat space limit, (u,v,w) be- 
comes (1, z 1 , z 2 ), where z 1 = x 1 + ix 2 and z 2 — x 3 + ix A 
parametrizing the R 4 space. In this limit, Eq.([T| can be 
written by 
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The R 2 version of Eq. the Laughlin wavefunction, 
is well known to be the unique ground state of the 
Trugman-Kivelson Hamiltonian [2]. To generalize, we 
shall show that Eq. ^ is the unique ground state of the 
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following Hamiltonian, 
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-xj), <5(fij) is the 



where fy = (x* - 

5-function, and Vy is the four-dimensional differential 
operator. To show Eq. ([2| to be the ground state of 
Eq. |3]), we first observe that for any pair i and j, Eq. ^ 
can be expressed by 



= [{z\ - z))F x + (z 2 - z*)F 2 f 
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where contain the symmetric part, including the 
even power of {z\' 2 — z 1 ^ 2 ). It is the property of the 
antisymmetry of the many-body wavefunction. Eq. ^ 
implies that the zero energy state should have (z\ 
z i) a ( z f ~ z ?) /3 where a + (3 = odd integers > 3. Eq. H 



can be easily checked to possess this property. Since 
Eq. ([3]) , describing the short-ranged repulsive interaction, 
is positive definite, Eq. ([2]), the Laughlin wavefunction 
in R 4 , is the ground state of Eq. ^jj. Moreover, using 
the similar argument in ref. [TJ |3] , the uniqueness of the 
ground state can be easily proved. We further note that 
the finiteness of the excitation gap of Eq. ^ is already 
shown in ref. E] , since the calculation was performed 
in the R 4 limit. Without elaborating too much technical 
detail, we remark that Eq. ^ is actually the flat space 
limit of the pseudo-potential Hamiltonian in CP 2 given 
in ref. pQ . 

Next, let us proceed to transform the theory to a lattice 
model. The dynamics of QHE in R 4 in the lowest Lan- 
dau level can be viewed as if there are two-independent 
non-commutative planes [H H]. One can choose them 
to be the x 1 -x 2 plane and the x 3 -x A plane. By applying 
the periodic boundary condition to the x 2 and x A direc- 
tions, the momentum in these two directions is quantized. 
Then, the single-particle wavefunction can be written by 
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where L 12 are the linear dimension in the x 1 ' 3 directions, 
m and n are the quantum numbers of the momentum in 
the x 2 and X4 directions respectively, Ib = y^hc/eB is the 

magnetic length, and a\ >2 — 77^-- These single-particle 
orbitals are Gaussian wave packages which localize in the 
x\ and 23 directions denoting by m and n. It defines a 
two-dimensional lattice with lattice constants ai, 2 , and 
m and n are the label of sites. Since L12 are finite, it 
automatically sets the upper bounds of m and n, that 
are L\ja\ and L 2 /a 2 respectively. Furthermore, 4> mn for 
different m and n are orthogonal. Therefore, 4>mn can be 
viewed as Wannier basis in this system. 



In order to restore the translational symmetry, periodic 
boundary condition has to be applied in the x 1 and a; 3 
directions. In this case, 4> mn is modified to be 
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Now, defining the field operator on this Wannier basis 
*0*0 = J2 mn c mn4>mn(x), where c mn is the annihilation 
operator of electron at the site (m,n), the second quan- 
tization of Eq. ([3| can be obtained easily as follows 
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f(l, I'; d, d)cl +l , tq+d ,<? p _ v q _ d ,c p+ i, q+ dCp-i, q -d (7) 
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f(l,l'-d,d') = 
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and Ki = di/ls- If the hopping ranges are measured 
in the unit of lattice constants, we found the hopping 
range ai/ Ki — Ib, which are the same in both directions. 
Furthermore, one can easily check that the hopping inte- 
gral f(l,l';d,d r ) in Eq. |7]) has a p4m group symmetry, 
known as the symmetry of the square lattice. It is certain 
that a\ = a 2 = a manifests the square lattice. Therefore, 
Eq. ([7]) describes a lattice model in the two-dimensional 
square lattice with periodic boundary. 

Eq. ([7]) is the two-dimensional generalization of the 
model proposed by Seidel et al. [5J. It is a pair hop- 
ping model with center-of-mass position conservation. In 
other words, it commutes with the following operators si- 
multaneously: 



U, 
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where i runs over all lattice sites which are labelled by two 
integers (x, y). Since the lattice has a periodic boundary, 
the system is translationally invariant. Namely, Eq. ^ 
also commutes with the translational operator T x and T y 
which translate the system by one lattice constant a in 
x and y directions respectively. However, Ui and Tj do 
not mutually commute. At filling 1/to 2 , they satisfy the 
following algebra 



T x U x Tl=U x e^, 



TyU y Tl 



(9) 
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where Li — m^a and U are integers chosen to be 
prime to to. Eq. ^ implies that the ground states 
cannot be simultaneously the eigenstates of T and U. 
One can choose the ground state to be the eigenstate 
of U satisfying U^e***, e^", E >= e*'^ 



and U y I e 



i<t>x e i<t> 



>,E > 



e**» , e^y , E 



,E > 
>. Then, 



TgTjjle 1 '^ , e l *» , E > are the orthogonal ground states for 
p and q to be the integer from to to — 1. In fact, the 
algebra in Eq. (|9| imply that the whole spectrum has at 
least TO 2 -fold degeneracy [7J. 

Now, let us analyze the ground state of Eq. Q of to = 
3. For k « 1, which can be achieved by taking Ib « L, 
the Wannier basis 4> mn highly overlap in space, since the 
lattice constant a decreases faster than Is- The ground 
state corresponds to the FQH states of Eq. ^ . Different 
from Eq. pj), the ground state of Eq. ^ is at least 9- 
fold degeneracy. The generate ground states are actually 
distinguished by the center-of-mass positions. As we shall 
show later, the ground state for small k indeed exhibits 
vanishing local order parameters. For n » 1, one can 
expand the Hamiltonian in terms of e~ K / 2 . We found 
Eq. can be well approximated by 
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where rij = cjcj and O(e~5- K e~5 K , e _3ft , e~2 K ) de- 
notes the positive energy for a pair gained by making the 
separation distance less than or equal to \/8a before or af- 
ter hopping. Because / > 0, Eq. ( |To| stabilizes a ground 
state of charge density wave (CDW) such that no pair is 
separated less than or equal to \/8a. We plot the ground 
states for N = 4 in Fig. ([T]). Given the center-of-mass 
position fixed, there are more ground state degeneracy in 
addition to the 9-fold degeneracy. For N — 2, the degen- 
eracy is 3, for N — 4 there are 5 shown in Fig. ([l^) to 
Fig. (Jljj), and for N — 6 there are 15, and therefore the 
total ground state degeneracy is 27, 45, and 45 respec- 
tively. In general, for N » 1, the degeneracy for fixed 



center-of-mass position is M = 3 



h-i 



1 with l\ 



and l 2 to be indivisible by 3, and the total ground state 
degeneracy is 9M. We stress that the result of the ad- 
ditional degeneracy is not an artifact of the truncated 
Hamiltonian. We performed the exact diagonalization 
calculation in the N = 2, 4, and 6 cases. We not only re- 
produce the analytical results given above but also obtain 
that the ground state degeneracy is independent of k. In 
the k » 1 limit, the uniform ground state can be ob- 



tained by the equally-weighted linear superposition of the 
degenerate ground states for a fixed center-of-mass posi- 
tion. The one in the N = 4 case is plotted in Fig. |]J). 
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FIG. 1. (Color online) (a) to (e) the ground states in the 
K » 1 limit in the N = 4 case. These five states have the 
same center-of-mass position, (f) the state of the equally- 
weighted linear superposition of the states from (a) to (e). 
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FIG. 2. (Color online) (a) the quasiparticle-quasihole pair can 
be created by shifting an electron up. The red curve arrow 
denotes the hopping of electron, (b) an example of the string 
excitation. The energy of the string excitation is independent 
of the length of the string, indicating deconfinement. The 
red ball denotes the location of the quasipaticle, and the blue 
circle is the one of the quasihole. 



From Eq.(10), it can be clearly seen that a excitation 



can be created if the distance between any two particles 
is equal to or smaller than y/8a. We can simply move an 
electron up by one site as shown in Fig. ([2^). Then, one 
can make a string excitation by moving more electrons 
up as shown in Fig. ^jp) . It turns that the energy of the 
excitations in Fig. (|2^) and Fig. ^p) is actually the same, 
because creating a string excitation does not make more 
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pairs to be separated less than or equal to >/8a. There- 
fore, the fractional excitations occur at both ends of the 
string. One end corresponds to the quasipaticle and the 
other end is the quasihole. The creation of quasiparticle 
and quasihole pair costs finite energy ~ e~ 2K . However, 
it does not cost energy to separate them. In other words, 
the quasiparticles and quasiholes are in the deconfined 
phase. 
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FIG. 3. (Color online) (a) a quaishole is a dislocation locating 
at (k, I), (b) a domain wall formed by the quasiholes 



The charge of the quasihole can be computed as the 
following. In a L x x L y lattice with open boundary condi- 
tion, let us create a quasihole at (k, I) shown in Fig. (|3^). 
It is nothing but inserting an empty site and create a 
dislocation in the system of CDW. Let us assume that 
its charge is e*. One can add more empty sites at row 
I so that a domain wall is created shown in Fig. ^jp). 
The charge of the domain wall is thus e*L x . Now, one 
can create one more domain wall by adding L x sites. If 
we add three domain walls, it is effectively equivalent 
to remove L x /3 electrons. Since three domain walls has 
charge 3e*L x , we obtain e* = e/9, where e is the elec- 
tric charge. For general m, the fractional excitation has 
e* = e/m 2 . It is easy to generalize to £>-dimensional 
integer lattice. Quasiparticle/quasihole is a dislocation 
behaving like a point-like particle in the D-dimensional 
lattice, where the fractional charge is given by e/m D . It 
is not surprising that the charge depends on the dimen- 
sionality, because the fractional excitation is a point-like 
particle. 

We believe that the FQH state at k << 1 can be adi- 
abatically connected to the charge-density-wave state at 
K » 1 without phase transition. In one dimension, Sci- 
del et al. has proven it true [6]. Here, we perform three 
calculations to support this conjecture. First, we checked 
that the ground state degeneracy is independent of n. We 
have performed the exact diagonalization to prove it true 
in the systems of 18 and 36 lattice sites. Second, we com- 



pute the following order parameter as function of k 

0(3^,0) = jy^e^cLcmn. (U) 

In the finite size system, Oi** ) is a function of L x and 
L y . In the inset of Fig. Q, we fix L x and show the 
scaling of Or 2* \ with L y . Then, taking the infinite L y 
limit, we plot the scaling of Q \ with L x . We found 
that the scaling of 0) with L x is almost identical to 
the ID case given by Seidel et al. [5]. The scaling shows 
very nice convergence. The smooth change of the order 
parameter indicates that there is no phase transition be- 
tween the k << 1 and k » 1 phases. Third, although it 
is not shown here, we also compute the energy gap for the 
system of 36 lattice sites. We found that the energy gap 
actually increase at around K — 0.5 where the order pa- 
rameter is about to deviate from zero observed by naked 
eyes and remains gapped in the whole region of crossover. 
We note that similar to the ID case the order parame- 
ter actually decreases exponentially below k < 0.5. With 
those three calculations mentioned above, we believe that 
the incompressible quantum liquid at k << 1 hosts the 
excitation of the charge e* = e/m 2 . 




FIG. 4. (Color online) In the inset, the scaling of with 
L y is calculated in the following cases L x x L v — 6 x 3, 6 x 6, 
and 6x9. The convergence is so fast that they are hardly 
distinguishable. Taking the value of O^jl i Q the infinite L y 
limit, the scaling of O-ijl with L x is computed for L x = 
6,9, 12, and 15. 

In summary, we show the existence of the fraction- 
ally charged excitation in a lattice hopping model in the 
square lattice. Our construction can be systematically 
generalized to higher dimensional square lattice, namely 
the integer lattice. The current theory answers positively 
that the fractionally charged excitation can in principle 
exist in any dimension. Our model is the first one to host 
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the ground state of incompressible quantum liquid in the 
integer lattice. 
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